The toric-boson model: Toward a topological quantum memory at finite temperature 
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We discuss the existence of stable topological quantum memory at finite temperature. At stake here is the 
fundamental question of whether it is in principle possible to store quantum information for macroscopic times 
without the intervention from the external world, that is, without error correction. We study the toric code in two 
dimensions with an additional bosonic field that couples to the defects, in the presence of a generic environment 
at finite temperature: the toric-boson model. Although the coupling constants for the bare model are not finite 
in the thermodynamic limit, the model has a finite spectrum. We show that, in the topological phase, there is 
a finite temperature below which open strings are confined and therefore the lifetime of the memory can be 
made arbitrarily (polynomially) long in system size. The interaction with the bosonic field yields a long range 
attractive force between the end points of open strings, but leaves closed strings and topological order intact. 
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I. INTRODUCTION 

In recent years, it has become clear that models of computa- 
tion based on quantum mechanics are radically different from 
classical ones,'-^ However, the quantum wave function of a 
macroscopic system rapidly decoheres because of the inter- 
action with the environment.^ Decoherence is responsible for 
the appearance of a classical world. "^-^ Therefore, if we want to 
process quantum information, we must preserve it from being 
spoiled by decoherence. 

In the past 15 years, a great amount of literature in quantum 
information has been devoted to the theory of error correc- 
tion, using several different paradigms. We encounter here 
a profound difference with respect to the classical case: en- 
coding classical information in a system can be robust, with- 
out the need for error correction (e.g., magnetic hard drives). 
On the other hand, the fragility of quantum states to decoher- 
ence seems to require quantum error correction. Is that really 
so? Can one build robust quantum memory at the physical 
level, just like one does with classical memory, or does quan- 
tum mechanics forbid it? 

For suitably symmetric models of the environment, there 
exist decoherence-free subspaces (DFSs)^i However, it is con- 
troversial whether these subspaces can exist in practice, be- 
cause the conditions to obtain DFSs are unstable.*^ Moreover, 
DFSs only exist for a specific model of noise, while we are in- 
terested in the question of principle whether a quantum system 
can be protected from a generic interaction with the environ- 
ment, without operating from the external world (i.e., error 
correction). The only restriction we shall assume about the 
way the environment interacts with the system is the locality 
of the interaction. 

If the only requirement is locality of the interaction be- 
tween the environment and the system, we can imagine en- 
coding information in some collective, topological degrees of 
freedom, and obtaining the desired robustness at the physical 
level. Physical systems that exhibit topological order seem 
to have the appropriate topological features^. The toric code 



is an example of a self-correcting quantum memory, that is, 
a quantum memory that is inherently robust against arbitrary 
(local) perturbations without need for error correction at zero 
temperature 

However, robustness against perturbations is not enough. 
What one needs is a system in which quantum information can 
be stored for macroscopically long times, in the presence of 
arbitrary local interactions with the environment at finite tem- 
perature T. There are not only errors due to virtual processes 
but, even more importantly, errors due to thermal excitations. 
For example, the toric code in two and three dimensions (2£'s 
and 3Z)s) is not stable i"-'^''^-''* [See Ref. 15 for a discussion 
on generic stabilizer codes.] The toric code in four dimensions 
is stablc j'^''^'^^ but interactions in four spatial dimensions are 
hardly realistic. 

One wonders whether quantum information is not only log- 
ically different from the classical one, but also physically dif- 
ferent in that it cannot be stored efficiently in a passive way 
(in less than four dimensions). Therefore, we are faced with a 
fundamental question: does a stable quantum memory exist at 
all or the laws of nature forbid it in three or less dimension? 

In this paper, we take a step toward realizing topological 
quantum memories in two spatial dimensions that are stable 
against decoherence towards an arbitrary environment. We 
obtain this result by adding to the toric code a local interaction 
between the defects and a bosonic field {toric-boson model). 
This result comes at the cost of a compensating defect energy 
whose strength diverges with growing system size. However, 
the protection arising in the toric-boson model is not trivially 
equivalent to removing thermal defects by introducing an in- 
finite energy cost by itself. In our model, the physical energy 
scales (namely, the excitation spectrum) remain finite, and it 
is the induced long-range interactions between defects that 
render the quantum state stable to thermal fluctuations. This 
is similar to what one usually does in high-energy physics, 
where infinite coupling constants are used (often needed) to 
obtain finite energy spectra. 

We show that the toric-boson model realizes a quantum 
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memory that can be stored for macroscopically long times, in 
the sense that the relaxation times scale (polynomially) with 
system size. While we fall short of realizing a stable topolog- 
ical quantum memory using local finite coupling terms in the 
Hamiltonian, our model shows that quantum relaxation times 
scaling with system size can be obtained in a system with a 
finite excitation spectrum, using local coupling terms only. 

II. TORIC CODE 

The toric code model"' is defined on a 2D L x L square 
lattice with periodic boundary conditions, and with spin- 1/2 
degrees of freedom on the = 2L^ edges. On every ver- 
tex s we can define the star operator Xs = Yljes '^^ 
every plaquette p the plaquette operator Zp = Yijep '^j' 
Hamiltonian of the toric code is 

HtoHc = -^(^X, + Y,Zp^ ■ (1) 

The operators Xg , Zp all commute and therefore thek eigen- 
values are good quantum numbers for labeling the states. Nev- 
ertheless they are not all independent because periodic bound- 
ary conditions imply Yl^ — Yip — ^- The two de- 
grees of freedom not labeled by the above quantum numbers 
are the two topological qubits defined by the algebra of the 
Pauli operators {Xl, ^ IljeL^ ^^j ' ^a/^ = UjeAh '^h 1 
{Xl^, Zmi} defined similarly, where Li, L2 and Mi, M2 are 
appropriately defined non-contractible loops around the torus. 

Two topological qubits can then be defined as the sub- 
spaces on which the two Pauli algebras [Zmi, Xl^) and 
Xl^) act, respectively. For example, one can choose 
the two basis states for the first qubit to be the eigenvectors 
|0) and |1) of Zmi, with eigenvalues +1 and —1. In this case, 
the operator flips the qubit {Xl.^ |0) = |1)), while the 
operator Zm2 can be used to change the phase of a qubit state 
[e.g.,ZM.(|0) + |l))/x/2-(|0>-|l))/x/2]. 

Since the model has a spectral gap A, the information en- 
coded in the topological qubits is stable against arbitrary lo- 
cal perturbations.'" However, it is fragile when interacting 
with an arbitrary environment at finite temperature (even if 
it is local) A heuristic argument for the fragility can 

be formulated as follows. Consider for concreteness Hi = 
g^f^iicTj (g) fj + cr| (8) fj) to be the interaction Hamilto- 
nian between the system and the environment. This inter- 
action will create pairs of defects (anyons) with probability 
P ^ exp(— A/ZcsT) and move them about in a Brownian 
motion. When two anyons travel around the non-contractible 
holes of the torus to recombine and annihilate, a logical error 
is produced, namely a non-contractible loop. [For example, 
a pair of Zp defects (generated by operators) that travel 
around a loop of type L2 will change the eigenvalue of Zmi , 
and VAfi.] Since the anyons propagate freely, the time scale 
for an isolated pair to induce an error is proportional to L^. 
However, the very same non-interacting nature that ensures 
free motion of the defects causes a non-vanishing density of 
them at any given temperature. Therefore, one can envisage 



faster collective processes whereby pairs stretch only as far as 
reaching anyons from neighboring pairs, and by annihilating 
them, they lead to longer overall propagation paths. The time 
scale to induce an error via these processes is roughly given 
by the time it takes for a pair to stretch over a distance equal 
to the average inter-particle separation, t^ci ^ exp(A/fcs7')- 
The relaxation time is therefore microscopic, because it does 
not scale with system size. 

The reason why end-points of strings (the anyons) can prop- 
agate freely is that here only the boundary of string-like ob- 
jects pays an energy, and a longer string does not have a larger 
boundary. So one way of confining these objects is to go to 
higher dimensions. In the 3D toric code, topological infor- 
mation is stored in non-contractible membranes and winding 
loops. Defects to these operators acquire the form of open 
strings and open membranes, respectively. While open strings 
still behave as in 2D, with end-points propagating freely, an 
open membrane costs energy that grows with the size of its 
boundary (a closed string), and it is therefore confined. Since 
only one of the two operators is preserved at finite temper- 
ature, the 3D toric code can only encode classical informa- 
jjQj^i3ji4 jj-^g quantum memory can be stored for times 

that are exponentially large in L, because both types of errors 
are membranes. "-'^-'^ 

One can think of confining the strings by adding a tension 
term to the Hamiltonian. There is a critical value of the cou- 
pling with this field such that the open strings are confined'^. 
Unfortunately, the tension term does not distinguish between 
open and closed strings. If open strings are confined, the sys- 
tem undergoes a quantum phase transition where topological 
order is destroyed, and with it the degenerate ground state en- 
coding the qubit is also lost. In fact, the phase in which the 
strings are confined can be mapped to the ferromagnetic phase 
of the 2D Ising model'**, and we can only encode classical in- 
formation in such system. 



III. TORIC-BOSON MODEL 

If we want to hinder the propagation of the anyons without 
affecting the loops and spoiling topological order, we must 
find an interaction that only couples the end-points of strings. 
A projector onto the subspace of open strings would allow us 
to give a tension only to open strings. But this term would be 
highly non-local, and therefore not realistic. Another possi- 
bility is to add a long-range force that makes the anyons at- 
tract, as it was suggested in Ref. [ll]. However, the problem 
is whether the needed attractive force can materiahze from a 
purely local Hamiltonian. 

As we show in the following, one can reach this goal with 
a local Hamiltonian by coupling the 2D toric code to bosonic 
fields, which will have the desired confining effect in the pres- 
ence of a completely generic bath acting on both systems. 
However, we shall achieve this at the cost of a compensating 
energy term in the Hamiltonian that scales with the size of the 
system. The bosonic fields will affect open strings by attract- 
ing the end-points, while they have no effect whatsoever on 
closed strings. In this way, they protect the topological qubit 
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while preserving topological order. 

Let Us and rip be the defect density operators [namely. 
Us — {1 — Xs)/2 and rip = (1 — Zp)/2, respectively] taking 
values and 1 on the stars / plaquettes of a square lattice (with 
spacing set to unity). This density will couple locally to the 
displacement field defined in the continuum. We consider the 
case of elastic waves (phonons) on the surface of the torus, 
two in plane and one out of plane. These elastic deformations 
are denoted by (y9(x, i), (/)(x, i) — ((/i^(x, t), 0*'(x, i)), where 
X takes values in the continuum [with centers of plaquettes 
labeled by Xp, and those of stars by x^]. 

The densities rig and Up couple to the out-of -plane phonons 
piezoelectric ally ((^)'^, and to the in-plane phonons in a ro- 
tationally invariant form (via dx<j>^ + dyip^). In the case of 
piezoelectric coupling, one has to take proper care of the cen- 
ter of mass (CM) mode. One can work with <y9k=o pinned to 
the origin (basically, the entire mass of the crystal behaves as a 
large classical object). This is physical, since what matters in 
the piezo coupling are the elastic deformations, and uniform 
translations should not affect the system; by working with co- 
ordinates relative to the CM, this is ensured. 

One can expand 93 (x, t) , (j)^ (x, t) using creation and annihi- 

lation operators '^k' ' ''k satisfying the harmonic os- 
cillator algebra as 



V3(x,t) = 

and 

</)Xx,t) = 




At) 



— ik-: 



(2) 




(3) 



with i — x,y, thus implying the canonical equal- time com- 
mutation relations 



[0^(x,t),0*(x',t)] = *<5(x-x'), 
[(p(x,i),93(x',t)] = i(5(x-x') 



(4) 



for the displacement fields. The Hamiltonian for the fields is 
given by 

i?boson = 5I'^k4«k+ (5) 

k/0 \^^0,i=x,y 

The interaction Hamiltonian reads as 

l—s, p 

(6) 

(notice that g^^ and have different dimensions), and the 
toric -boson Hamiltonian can thus be written as 



H = Ht, 



(7) 



After substituting Eqs. (O and ^ into the interaction 
Hamiltonian (|6]l, and introducing new creation and annihila- 
tion operators 



"k 


= «k 




= ^k 



9uj pH'k) 
ign^i P^(k) 



(8) 

{^ = x,y), (9) 



the total Hamiltonian H of the system can be written as 



(Of flW 



H = i/toric + ^ c^k + ^k/3r^A 

kT^O k^O,i=x,y 



9Z 



2^1 



(10) 



where p(k) = J^e nie^^'^'^ is the Fourier transform of p{x) = 
rig (5(x — Xf). Here we have assumed that the defects are 
smeared over the corresponding plaquette / star, so as to regu- 
larize (5(x) and remove ultraviolet divergences. 

Notice that the first two terms in Eq. ( fTOl ) lead to precisely 
the same energy levels as those of the bosonic modes in the 
absence of the defects, [Eq. (|5]l]. This holds true even for ther- 
mal averages in mixed states. Therefore, the phononic fields 
added to the toric code Hamiltonian give rise to an effective 
interaction potential Vd given by the last term in Eq. ( fTOl ), ir- 
respective of temperature. Notice that Vd is always negative, 
independent of the arrangement of the defects, and thus it fa- 
vors defect proliferation, which will need to be appropriately 
counter-balanced later on. 

Consider the case of 7V^ = Mx + ■N'z defect pairs, at po- 
sitions Xa, with t — 1, . . . , 2Mx labeling the sites with star- 
type defects and I — 2Mx + 1, ■ • ■ , 2M labeling the sites with 
plaquette-type defects. The effective interaction potential Vd 
can then be written as 



Vd (Xi, . . . ,X2AA) = - 

I 


2A/' ^ 

,t' = l kT^O 


X 


' 9I 




[2^1 
1H 


2^1 \ 



(11) 



K;(xf - X£/) + Vh(xc - X£() I] 



£,£' = 1 



where 



K;(X£ - X£' 

Vn{xi - X£/ 



5c£_gjk-(xj-Xj/) 



-E 

1 v^g^lkp 

v^ ^l' 

k/0 k 



ik-(x^— x^/) 



(13) 



(14) 



For acoustic phonons with dispersion o^k — vj^\, and 
i^k = Wf2|k|, the potential V^iv) becomes Coulombic in d- 
dimensions. 



9l 1 - 



k^^O 



(15) 
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Namely, the potential generated by a point-like charge in a 
uniform compensating background distribution is 



Vu{r) 



gt T{d/2) 1 I |2-d 
^7" 27r''/2 2-d I' I 



% 1 Inlrl +C 



C d^2 



d= 2, 



where the uniform compensating background is irrelevant due 
to the absence of the k = term in Eq. ST3[ . For the 
same reason, JV^{x.) d'^x — 0, which fixes the constant 
C = —g'^/{2TTv'^) In^i, with f/, of the order of the system 
size L in two dimensions (the actual value will depend on de- 
tails about the shape of the system). 

For > 0, the potential K;(r) acts therefore as 
as a "gravitational" potential with Newton constant G = 
{2v^'K'^^^)^^r {d/2), attracting particles of mass g^^. 

The in plane phonon contribution amounts to 



Vr 



<5W(r)--^ 



and we can finally combine the two contributions into 



1 

'2^ 



In- 



V' 



2-E 



(16) 



(17a) 



(17b) 



(17c) 



where we separated for convenience ln(|r|/^i) into 
ln(|r|/a) — ln(^i/a) (a being the lattice constant, |r|/a > 
1, Vr). The first term is the desired confining interaction 
(which we discuss below), while terms (I17bb and ( I17cb are 
by-products that need to be appropriately taken care of. By 
tuning the ratio between gu/v^^ and gn/vn one can obtain 
that the two parts of Eq. (I17bb cancel each other out exactly, 
at the cost however of generating a system-size scaling cou- 
pling constant in Eq. ( I17cb : 



9R = 9l}^^^^ 
Uo v^, 2tt a 



InL. 



(18) 



Therefore, the phonon-mediated effective interaction between 
the defects can be written as 



Vd (Xi, . . . ,X2A/-) 

9l 1 



2 , 2M 

v^, 47r ' a 



niUi, In 



2 1 / t \ 



e,i'=i 



2 1 2AA 

= — — } mnt In- 

u,2 47r ^ 



-4— fyin^^2AA,5(2)(0), 
47r V a 



(20) 



where the (5-function has been regularized at short distances 
just like one does for the Coulomb potential. 

The remaining unwanted term - coming from Eq. ( I17cl ) 
- plays the role of a chemical potential that favors thermal 
defects, and whose strength scales with the size of the system. 
In order to preserve the fact that the defect density vanishes 
at zero temperature, we must neutralize this contribution by 
adding a compensating energy cost for the defects that also 
grows with system size. 

Since this unwanted term scales linearly in the number of 
defects M, it can be removed by the further addition of a 
chemical potential to the Hamiltonian of the system. Indeed, 
the role of the in plane phonons in our model is to prevent 
the arising of undesired defect-favoring long-range interac- 
tions that scale with the square of the number of defects M"^, 
and could not be removed by the addition of a local one-body 
term in the Hamiltonian. 

In conclusion, we obtain a purely gravitational defect inter- 
action of the form 



gl 1 



2A/' 



(xi, . . . ,X2^) = > ning.ln 

u, , 47r ^ — ' a 

" i,i'=i 



(21) 



An important comment is in order While we did add an 
infinite I? In L) coupling term to the Hamiltonian in or- 
der to arrive at the effective interaction shown above, this is 
not equivalent to trivially introducing an infinite energy cost 
for the defects, which would forbid their presence altogether. 
In fact, the energy to add a pair of defects in our model re- 
mains ^n;fe (the environment acts locally on the system and 
therefore defects are always produced close to each other, and 
move one step at a time). 

The situation is similar in high energy, where the bare 
Hamiltonian often has diverging coupling constants, but it car- 
ries nonetheless physical meaning as the spectral gaps are fi- 
nite. 



A. Robustness of quantum memory 

Once we obtain an effective long-range "gravitational" at- 
traction between the defects, we can argue for the stability of 
the quantum memory. It is well known that the partition func- 
tion of M particles of mass m interacting via a 2D gravita- 
tional potential u(r) — Grr? ln(r) becomes divergent below 
a finite temperature = MGm? /A (see e.g., Ref. 20). The 
gravitational forces overcome the entropic contribution and 
lead to a collapse, where all the particles coalesce to a single 
point. In the toric-boson model, where the massive particles 
correspond to the defective stars and plaquettes, their number 
is controlled by the chemical potential due to the ratio A/T. 
Therefore, given that the number of particles is topologically 
constrained to be even, the equilibrium density of defects van- 
ishes below the finite temperature T^,{J\f = 2) — Gw? /2. 

In this temperature regime, thermal processes connecting 
two different topological sectors are essentially limited to cre- 
ation, diffusion and annihilation of individual pairs of parti- 
cles. Such processes must overcome the energy barrier of tak- 
ing one particle around the whole system, which scales with 
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the logarithm of the system size, thus yielding a characteris- 
tic macroscopic time scale tfoi ~ exp[aln(L)/r] ~ i"/-'" 
which is polynomial in system size L [poly(L)]. 

IV. ROBUSTNESS OF THE TORIC-BOSON MODEL 

An essential ingredient to obtain the effective gravitational 
potential between the defects lies in the gaplessness of the 
phonon modes involved in the system. In order to ensure ro- 
bustness against thermal fluctuations in the toric -boson model, 
one should therefore also ensure that the action of the environ- 
ment on the phonon modes does not open a gap (e.g., via in- 
troducing effective interactions between the different bosonic 
modes). 

In our case, the gaplessness (wk ^ |k|) of phononic ex- 
citations is protected by translational symmetry. Even if the 
elastic medium is inhomogeneous, it satisfies the equation 
of motion m(x) 0(x, t) = J d'^y K{x, y) 0(y, t), where the 
condition J df^y i^(x, y) = follows from uniform trans- 
lations 0(x, i) = const.; this condition on the i4r(x, y) is 
similar to that in a master equation for stochastic systems^', 
with the difference that the wave equation is second order. 
It thus follows from generic arguments on diffusion in in- 
homogeneous media that the dispersion of the phonons for 
long wavelengths must be linear, ^ |k|. Notice that this 
conclusion holds even if the inhomogeneous elastic medium 
breaks translation symmetry [in the form of an inhomoge- 
neous m(x) and _fir(x, y)]. The key point is that, at long 
distances, acoustic sound waves propagate in the effective 
medium. Notice also that to attempt to pin the displacement 
fields with a rigid substrate potential is an illusion, as one must 
consider the elastic properties of the substrate itself. The lin- 
early dispersing acoustic modes in the effective medium at 
large distances are responsible for the the long-range behavior 
of w(r), which remains unchanged after tracing out the envi- 
ronment. [One can bypass the issue of the linear dispersion 
and compute directly Va, (r) in the inhomogeneous medium: 
Vh (x' , x) = dt P\Y (x' , x; t), where Pw (x' , x; t) is the 
probability to diffuse from x to x' in time t, computed using 
the transition rate W{x, y) — K{x, y) /m(x).] 

V. CONCLUSIONS AND OUTLOOK 

In this paper, we have shown how a system with topologi- 
cal order like the toric code can be made robust to local inter- 
actions with a generic environment below a finite temperature 
T*, if the defects in the system are coupled to a phononic field. 
The defects feel a long-range attractive force that prevents 
them from moving freely around the torus. The long-range 
force is not spoiled by the interaction of the bosons with the 
environment as long as the dispersion relation is preserved. 
In our model, phonons are protected by translational symme- 



try and the dispersion relation yields a logarithmic potential 
resulting in a relaxation time scaling as poly(L). Therefore, 
quantum information can be stored for arbitrarily long times 
scaling with the size of the system. 

In the model described here, we need to introduce compen- 
sating terms that scale as In L, L being the linear size of the 
system, and thus diverge in the thermodynamic limit. How- 
ever, we do so while preserving the finite nature of the low- 
energy spectrum of the system. We believe that this is a step 
farther with respect to the simple addition of an infinite en- 
ergy cost for the defects in the Hamiltonian. Our model shows 
that, in principle, quantum mechanics allows for the preserva- 
tion of quantum information for arbitrarily long times, scaling 
with the size of the system. 

There are a number of open problems left to investigate. 
First and foremost, can one construct a model with similar 
physics where all the interactions are not only local but also 
finite! 

Second, for reliable quantum topological memory, poly(L) 
is not ideal. It takes in general polynomial time to operate on 
a topological qubit (e.g., the read-out at finite temperature''*). 
Our system would be viable only in the low temperature limit, 
where the polynomial relaxation time scale ^rei ~ I/"^^ is 
much larger than the (temperature independent) polynomial 
operational time scale. A stronger quantum memory could 
be realized in a system that exhibits exponential relaxation 
time scales t^ei ~ exp(L). This could be obtained using a 
bosonic field with protected gapless excitations with disper- 
sion uj ~ There are other topologically ordered systems 
with the right cubic dispersion relation, whose gaplessness is 
protected by topological order itself, resulting in a very strong 
protection^-. Therefore, a toric code interacting with the ap- 
propriate system would feature exponentially large relaxation 
times. We think that the robustness of quantum memory in 
such system is the true meaning of the full survival of quan- 
tum topological order at finite temperature. 



VI. ACKNOWLEDGMENTS 

We acknowledge important discussions with M. Hastings, 
X.-G. Wen and D. Lidar. We are particularly grateful to 
M. Hastings for pointing out some of the undesired terms that 
arise in our Hamiltonian. Such terms would have led to phase 
separation effects in the system, and we had missed them in 
a previous version of the manuscript. Research at Perimeter 
Institute for Theoretical Physics is supported in part by the 
Government of Canada through NSERC and by the Province 
of Ontario through MRJ. This work was supported in part by 
EPSRC under Grant No. GR/R837 12/01 (C. C), and by a 
grant from the Foundational Questions Institute (fqxi.org), a 
grant from xQIT at MIT (A. H.). A. H. and C. C. thank the 
BU Condensed Matter Theory visitors program for its hospi- 
tality. 



' D. Deutsch, Proc. R. Soc. London, Ser A400, 97 (1985); D. P. Di- Mod. Phys. 68, 733 (1996). 

Vincenzo, Science 270, 255 (1995); A. Ekert, and R. Josza, Rev. 



6 



^ M. Nielsen, and I. L. Chuang, Quantum Computation and Quan- 
tum Information, (Cambridge University Press, Cambridge, UK 
2000). 

^ I. L. Chuang, R. Laflamme, P. W. Shor, W. H. Zurek, and Science 

270, 1633 (1995). 
* W. H. Zurek, Phys. Today 44, No. 10, 36 (1991); W. H. Zurek, 

S. Habib, and J.P Paz, Phys. Rev. Lett. 70, 1187 (1993); 

W. Zurek, Phys. Rev. D 26, (1982). 
^ D. Giulini, E. Joos, C. Kiefer, J. Kupsch, I. O. Stamatescu, and 

H. D. Zeh, Decoherence and the Appearance of the Classical 

World (Springer- Verlag, Berlin, 1996). 
^ P W. Shor, Phys. Rev. A 52, 2493 (1995); A. Steane, Phys. Rev. 

Lett. 77, 793 (1996); E. Knill, R. Laflamme, and W. H. Zurek, 

Proc. R. Soc. London, Ser. A 454, 365 (1998); D. Gottesman, 

PhD thesis , California Institute of Technology, 1997. 
^ P Zanardi, and M. Rasetti, Phys. Rev. Lett. 79, 3306 (1997); 

D.A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. Lett. 81, 

2594 (1998). 
^ R. Ahcki, Physica A 150, 455 (1988). 

' X.-G. Wen, Quantum Field Theory of Many-Body Systems, (Ox- 
ford Univ. Press, Oxford, 2004); Phys. Rev. B 40, 7387 (1989); 
X.-G. Wen, and Q. Niu, Phys. Rev. B 41, 9377 (1990); Int. J. 
Mod. Phys. B 4, 239 (1990); Adv. Phys. 44, 405 (1995). 

" A. Kitaev, Ann. Phys. 303, 2 (2003). 



" E. Dennis, A. Kitaev, A. Landahl, and J. Preskill, J. Math. Phys. 
43, 4452 (2002). 

Z. Nussinov, and G. Ortiz, Phys. Rev. B 76, 174416 (2007), and 

arXiv:cond-mat/0702377 (unpublished) 
" C. Castelnovo and C. Chamon, Phys. Rev. B 76, 184442 (2007); 

Phys. Rev. B 78, 155120 (2008). 
'* R. Alicki, M. Fannes, and M. Horodecki, J. Phys. A: Math. Theor. 

40, 6451 (2007); J. Phys. A: Math. Theor. 42, 065303 (2009) 

A. Kay, and R. Colbeck, arXiv:08 10.3557 (unpubhshed); 

S. Bravyi, B. Terhal, New J. Phys. 11, 043029 (2009). 

R. Alicki, M. Horodecki, P. Horodecki, R. Horodecki, 

arXiv:08 11.0033 (unpubhshed). 
" S. Trebst, P Werner, M. Troyer, K. Shtengel, and C. Nayak, Phys. 

Rev. Lett. 98, 070602 (2007). A. Hamma, W. Zhang, S. Haas, and 

D. A. Lidar, Phys. Rev. B 77, 1551 II (2008); 

A. Hamma, and D. Lidar, Phys. Rev. Lett. 100, 030502 (2008). 
" G. D. Mahan, Many-Particle Physics, 2nd ed. (Plenum, New York, 

1990, Sec. 1.3 D). 

A. M. Salzberg, J. Math. Phys. 6, 158 (1965); P -H. Chavanis, 
Eur. Phys. J. B 57, 391 (2007). 

T. Nakayama, K. Yakubo, and R. L. Orbach Rev. Mod. Phys. 66, 
381 (1994) 

X.-G. Wen, private communication. 



